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1. INTRODUCTION 

Smart grids are a set of technologies, concepts and approaches, allowing the integration the 
generation, transmission, distribution and use into one internet by full use of advanced sensor measurement 
technology, communications technology, information technology, computer technology, control technology, 
new energy technologies [1]. However, Smart Grid uses digital technology to control grid and choosing the 
best mode of power distribution to reduce energy consumption, reduce costs, increase reliability and also 
increase transparency in the network. Therefore, the system intelligent will have will have a significant 
impact in the fields of finance and economics of the power industry [2]. Although, the traditional network is a 
one-way network in which the electrical energy produced in power plants is channeled to consumers without 
information to create an automated and distributed network of advanced power supplies. 

The unit commitment problem plays a significant role in optimizing the cost of generating electrical 
power by planning production units based on the allocation of the production cost of each unit and the actual 
output power [3]. They involves scheduling the on/off states of generating units to minimize the operating 
cost for a given time horizon. The committed units must meet the systems fore-casted demand and spinning 
reserve requirement at minimum operating cost, subject to a large set of operating constraints. 

The UC problem, one of the most important tasks in short-term operation planning of modern power 
systems, has a significant influence on the secure and economic operation of power systems [4]. 
Optimal commitment scheduling cannot only save millions of dollars for power companies; it also ensures 
system reliability by maintaining the proper spinning reserve. 
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ED is also applied in the integrated system for scheduling power plants. A few methods have been 
published to solve the ED problem and Optimal Power Flow (OPF). Researchers have published a few 
methods to solve ED and OPF problems. Direct method is accurate and very simple but limited by the 
quadratic objective function [5]. 

The economic dispatch (ED) is one of the power management tools that are used to determine real 
power output of thermal generating units to meet required load demand. The ED results in minimum fuel 
generation cost, minimum transmission power loss while satisfying all units, as well as system 
constraints [6], [7]. 

The rise of energy demand and insufficient of energy resources are required for quality and secured 
dispatch [8]. A well-coordinated and optimized power system operation helps in satisfying Economic 
Dispatch (ED) among users of power networks. Hence, studies need to be conducted in order to analyze and 
develop new tools so that the optimization issues in ED could be overcome. Basically, the principal objective 
of load dispatch is to minimize the total fuel cost while satisfying the requirements of some important 
operational parameters. In today’s environment, efficient load dispatch requires not only to schedule the 
power generation at the least cost but also to consider other performance factors to be optimized in power 
flow over the networks. The obligation of social attentions has influenced the reduction of energy 
conservation and pollution emission produced by power plants [9]. 

Facing the electric bulimia experienced by the world and as an urgent and efficient solution is 
sought, it is essential to optimize the cost of producing electricity. As such, tiny cuts costs conceal huge 
potential savings, this is part of this paper, and we look at the overall optimization purposes known economic 
load distribution (OPF) or economic dispatch (ED) [10]. 

The ED is a static problem is to say we must define at a given powers generated by each power plant 
to power a load as economically as possible. To solve this problem the optimization methods are used. 

Conventional optimization techniques [11], [12]. Have long been applied to solve the ED problem 
such as Quadratic Programming [13], [14]. Linear programming [15] sequential approach with a matrix 
framework (SAMF) [16]. Modified Lambda-iteration method [17], Newton Raphson and Lagrangian 
multiplier (LM) algorithms [18], in the real-design cases, the number of decision variables (i.e. power units) 
of the ED area are very large. The objective criterion to be minimized could also have too many local 
minimum which might not lead to the minimum cost and the best generation schedule of power system units. 
Therefore, efficient search algorithms are needed. 

Nature-inspired metaheuristic search algorithms gain popularity due to their promising performance 
on solving many real-world optimization problems which are complex, nonlinear and multi-model. In the 
past two decades, the literature of metaheuristic search has expanded extensively. 

Some of the well-known metaheuristic approaches are Genetic Algorithms [19], Genetic 
Programming [20], [22], Particle Swarm Optimization [23], [24], Simulated Annealing [25], Artificial Bee 
Colony (ABC) [26], Cuckoo Search [27-28], crow search algorithm [29]. 

The rest of this paper is organized as follows. Section 2 contains the problem formulation of the ED 
and UCP. Section 3 briefly presents the basics of ES and CSA. Section 4 proposes the binary eagle strategy 
based crow search algorithm to solve Economic dispatch and UCP. Section 5 provides the computational 
results. Finally, Section 6 outlines the conclusions. 


2. PROBLEM FORMULATION 

The scheduling problem of generators solved ideally by acquiring exhaustive trial of all solutions 
and best solution is chosen amongst them. All possible units supplying a load and reserve requirements 
would be tested and choose the optimal solution that have the minimum operating cost [30] The generating 
units’ output power with system constraints over a time period T and startup/shut down times at each step 
required to scheduling problem of generator. The running cost significant term of a thermal unit is the output 
power of the committed units [31] the fuel cost, FCi is represented in a quadratic form of output power in a 
time interval given in (1) 


Fp = a5 F,(P,) = ee aj + b;P, + c,P? (1) 


Where a, b, c are cost coefficients of unit and p is the unit generating power. The start-up cost (SC) 
calculation depends on the treatment strategy for a thermal unit during down time periods and an 
exponentialcost curve shown in (2) is its representation. where o, 4, tT is the hot startup cost, the cold startup 
cost and the cooling time unit constant and T, is the time at which the unit has been turned off so The total 
production cost, F is the sum of the operating, startup and shut down costs for all the units illustrated in (3) 
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fi f ‘i Tofi 
SCi =oi + 6i * {1 — exp (-=) (2) 
Fr = Yi_, UN, FCi,t + SCi,t + SDi,t (3) 


Where N is the number of generators and different load demands number is T at estimated 
commitment, SD is the shutdown cost. 
Some constraints should be taken into consideration to minimize F as: 
(i) Power balance equation is given by (4): 


dik Pi Ui = (Pp + PF) (4) 


Where PD is the load demand and PL is the power loss of the system. 
(ii) The hourly spinning reserve (R) is given by (5): 


x=N, Pi,max Ui = (Py +P.) +R (5) 
(iii) Unit rated minimum and maximum capacities as in (6): 
Pumin s PB s Pimax (6) 


The initial conditions of each unit and minimum up /down (MUT/MDT) time limits of units are 
given by (7) and (8) respectively. 


(Ti — MUTI) * (Uy-1,1 — Uri) 2 0 (7) 
Cie y — MUTI) * (Uy-1,1 — Uri) 2 0 (8) 


Where the unit off/on time is Toff/Ton the and the unit off /on [0, 1] status is Ut,i. the enhancement 
of ELD problem is represented by (9): 


Fr = Dik F(P) DE aj + bP, + GP? (9) 
Subject to the equality and inequality contraints are given by (10) and (11) respectively. 
M1 Pi = (Pp +P.) (10) 


Piumin < B < Pimax (11) 


3. OVERVIEW OF EAGLE STRATEGY AND CROW SEARCH ALGORITHM 
3.1. Egale Strategy 

Eagle strategy is a two-stage optimization strategy was presented by [32]. This algorithm mimics 
behavior of eagles in nature. In fact, eagles use two different components to search for their prey. The first 
one is a random search performed by flying freely and the second one is an intensive search to catch prey 
when they see them. In this two-stage strategy, the first stage explores the search space globally by using a 
Levy flight: if it finds a promising solution, then an intensive local search is employed using more efficient 
local optimizer, such as hill-climbing and the down-hill simplex method. Then, the two-stage process 
commences another time with new global exploration, followed by local search in a new area. One of the 
remarkable advantages of such a combination is to use a parallel balance between global search (which is 
generally slow) and a rapid local search. There i s another advantage that is called a methodology or strategy, 
not an algorithm. In fact, there are different algorithms that can be used at different times and stages during 
iterations. The main steps of the ES are outlined in Algorithm 1. 
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Algorithm 1 Eagle strategy 

Objective function f(x) 

Initialization and random initial guess x‘~ 

While (stop criterion ) do 

Global exploration by randomization (e;g ;levy flights ) 
Evaluate the objectives and find a promising solution 
Intensive local search via an efficient local optimizer 
If (a better solution is found ) then 

Update the current best 

End if 

Update t=t+1 

End while 





mere OAAYNANANMNHWNK 


SD Siete ee 





3.2. Crow Search Algorithm 

The crow search algorithm (CSA) is a new population-based stochastic search algorithm recently 
proposed by [33]. The CSA is a newly developed optimization technique to solve complex engineering 
optimization problems [34-35]. It is inspired by the intelligent behavior of crows.The principles ofCSA are 
listed as follows [33]: 

Crows live in the form of the flock. 

Crows memorize the position of their hiding places. 

Crows follow each other to commit thievery. 

Crows protect their caches from being pilfered through probability. 

Following the above assumptions, the core mechanism of the CSA consists of three basic phases, 
namely initialization; generate a new position, and updating the memory of crows. At first, the initial 
population of crows represented by n dimension is randomly generated. At iteration t, the position of crow is 
specified by xit = Ese x a hes xy" land it is assumed that this crow has memorized its best experience thus 
far in its memory m!t = [mi me; are m}'|To generate a new position, crow i select randomly a crow j, 
for example, from the population and attempts to follow it to find the position of its hiding place (mj ). In this 
case, according to a parameter named awareness probability (AP), two states may happen: 

State 1: Crow j does not know that crow i is following it. As a result, the crow i will determine the hiding 

place of crow j. 

State 2: Crow j knows that crow j is following it. As a result, to protect its cache from being pilfered, 

the crow j will fool crow i by going to another position whitin the search space. 

According to States | and 2, the position of the crows is updated as follows: 
Stet ei xiiter+1 +1; X fliiter x (miiter _ xbiter) (12) 


12 Ap?iter A random positonofsearch spaceotherwise 


Where rj is a uniformly distributed fuzzy number from [0; 1] and AP!!*eT denotes the awareness probability 
of crow j at iteration iter. Finally, the crows update their memory as follows: 


miitert1 — ail if f(xilter )is better than f(miiter ) Pe 
miter otherwise 
Where f( - ) denotes the objective function value. It is seen that if the fitness function value of the new 


position of a crow is better than the fitness function value of the memorized position, the crow updates its 
memory by the new position. The above process is repeated until a given termination criterion (itermax) is 
met. Finally, the best solution of the memories is returned as the optimal solution found by the CSA. 
The main steps of the CSA are outlined in Algorithm 2: 
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Algorithm 2 Crow Search Algorithm 





di Randomly initialize the position of a flock of (Np ) crows in the search space. 
2 Evaluate the position of the Crows 

3: Initialize the memory of each Crow 

4: While (iter< iternax) do 

5: for i= 1: to Npdo 

65 Randomly choose one of the crows to follow (for example, j) 
Ts Define an awareness probability 

8: if (r,=AP/”) then 

9: xiitert] iter sy. ofl x(x ‘iter ) 

10: else 

ll: x'+l—a random position of search space. 

12: end if 

13: end for 

14: Check the feasibility of new positions 

15:; Evaluate the new position of the Crows 

16: Update the memory of crows 

Li end while 





4. BINARY EAGLE STRATEGY BASED CROW SEARCH ALGORITHM FOR ED AND UC 

The binary ES-CSA is used to optimise the unit-scheduling problem in the first step, and the 
Lambda-iteration method [36] is used to solve the economic load dispatch problem in the second step. 
These two steps run iteratively until the algorithm meets the stopping criterion. Optimising the first sub 
problem of unit-scheduling is more difficult than the other sub-problem of ELD. So this paper mainly 
discusses how to model BESCSA for the first sub-problem, and the second sub-problem is solved by the 
traditional Lambda-iteration method. These two sub-problems are optimised iteratively until the algorithm 
meets the stopping criterion. The (13) are transfer from continues to binary space using the following 
equations : 


xiter = {1 ifs(eile" ) 2 randO, (14) 
0 if otherwise 


Where s(x@ter) = = y=1 +e%" and rand() is a random number from uniform distribution [0; 1] and 


xtiter is the updated binary position at iter iteration. 


4.1. Solution Representation and Initialization 

Before using the proposed binary ES-CSA to solve UCP, the representation of a crow must be 
defined. A crow is also called an individual. Hence, we defined each unit on/off (or 1/0) status as a gene, 
all available unit status at each hour make up a sub-chromosome, and there are H sub-chromosomes over the 
time horizon H comprising an individual. An individual would display the unit commitment schedule over 
the time horizon H. The on/off schedule of the units is stored as an integer-matrix U with dimension N G H. 
A matrix representation of an individual in the population is shown as follows: 


ut uz uz ... ul 

us us us sia Eee 

U = 7 ~ a a 
2 ; 

uy urs we ... ul 


Where uhi is unit on/off status of unit i at time h (uhi = 1=0 for on/off). 

In the initialization process, a set of individuals is created at random. For the complete N P 
population, the candidate solution of each individual Uj; (j = 1; 2; :::; N P ) is randomly initialized. 
The position uhi of each crow Uj is generated using a uniform distributed random function, which generates 
either 0 or | and they are equally likely. 


4.2. Generate New Solutions 

As mentioned above, the ES is a two-stage strategy, and we can use different algorithms at different 
stages. In the first stage, ES uses the so-called Levy flights, which represent a kind of non-Gaussian 
stochastic process whose step sizes are distributed based on a Levy stable distribution to generate new 
solutions. When a new solution is produced, the following Levy flight is applied: 
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xhiter+1 = ybiter + @ @ Levy(A) (15) 


Here, is the step size that is relevant to the scales of the problem. The product means entry-wise 
multiplications. Levy flights essentially provide a random walk while their random steps are drawn from a 
Levy distribution for large steps: 


Levy(a) =u = t*,1<A<3 (16) 


In this paper, we will use the Mantegna algorithm [37], which is one of the most efficient algorithms used to 
implement Levy flights. We assume that Levy (A ) = s, so the formula can also be described as follows: 
By using Mantegna’s algorithm [38], the step length s can be calculated as follows: 


xliter+1 — yiiter 4 q Os (17) 


By using Mantegna’s algorithm [26], the step length s can be calculated as follows: 


—H1 
s=i /p (18) 


lvl 


Where wand v draw from the normal ditrutions respectively. that is: up ~ N (0,02),v ~ N (0,02), and 
r(a+p).sin@® 2 ; ; 
= re oy, = 1 Here 0< B < 2 and (:) is theGamma function. 
r(&*) 2° 
For the second stage, we can use the crow search algorithm (CSA) for the intensive local search. We know 
the CSA is a global search algorithm, but it can easily be tuned to do an efficient local search by limiting new 
solutions locally around the most promising region. As mentioned above, in the CSA, there are two specific 
parameters: awareness probability (AP ) and flight length (f 1). Small values of AP intensify the local search, 
while large values result in a global search. Hence, the CSA can easily be used as a local optimizer by setting 
the awareness probability to very small values, and for good performance, we choose the flight length f 1 = 2. 
Such a combination may produce better results than those using pure CSA. 

In UCP, binary numbers 0 and | are used to indicate the unit status (i.e., OFF or ON). The proposed 
startegy is essentially a real-coded algorithm, and therefore some modifications are needed to enable it to 
deal with the binary variable (i.e., 0 and 1) optimization problem. 


Oy O,, are calculated as follows ( 


5. RESULTS AND DISCUSSION 

In this section, we present the results obtained based on ES-CSA for solving the economic dispatch 
problem and unit commitment problem and compare this results with the CM (Conventional Method) [17] 
and Particle Swarm Optimization [39]. A 10 unit’s power unit system to explore our idea on using ES-CSA 
to find the optimal set of power generation of the system. ES-CSA will be used in this paper to solve the 
problem of economic dispatch and unit commitment. The programs are developed in MATLAB 7.9 
environment. 

The adopted system is expected to produce demand power of 1500 MW. The tuning parameters for 
ESCSA are given in Table 1; the Table 2 shows the cost coefficient of the 10 generators, under study, 
while the matrix is the loss coefficient matrix of the 1Ounits power system. the generating unit data of the test 
system are given in Table 3 From the results of Table 4, we notice that ES-CSA give us the same production 
cost, and CM gives a slightly lower cost of $ 0.7 / h, ES-CSA gives us a good production cost and good 
accuracy. In the meantime, we examine the variation in the total fuel cost of test system with evolutionary 
generation numbers. For differenttest systems, the convergence processes of the best solution in the30 trials 
are listed in Figures 1 and 2. From Figure 1, it is easy to see the ES-CSA has satisfactory con-vergence and 
the algorithm escaped from the local optima at thelater iterations. It proved that the stochastic searching 
mechanism of ES-CSA, which is conducted by gravitational forces among agents, is efficient. And the 
proposed mutation strategies improved theperformance of ES-CSA. 

In Figure 2, we show the convergence of the metaheuristic search process based on ES-CSA in both 
the best and average cases. To see the difference between our new approach and another known method, 
we will compare the production cost found by ES-CSA to that found by PSO [40]. 
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In this case, we will test the operation of ES-CSA. For this, we will use a simple network of 26 
nodes with 10 production units. The total demand of the network is equal to 1500 MW and loss coefficients 
are as follows: 








49 14 15) 15 £6 1:7. 7 2S 19 20] 
14 45 16 16 17 15 15 16 L8 18 
15 1:6 39. 10 12. 1.2: 14 24 £6 1.6 
15 16 10 40 14 1.0 1.1 12 14. 215 
B ~ 10-5 16 17 12 1A 3S 1.1 13 13 15 16 
ae Ly tS 12 1.0 11 336 12> 12 TA: ES 
17 15 14 LT 13 12 38 1.6 bs 1.8 
Ls te 14 12 13 £2. 16 40 TS 1.6 
i913 16 14°15 14 16 75 42 19 
{2.0 18 16 15 16 15 18 16 19 44] 


The simulation results are presented in Table 4. 


Table 1. Parameters of CSA AND ES-CSA[41]. 








Algorithms/parameters AP fl B 
CSA 0.2 2 - 
ES-CSA 0.2 2 15) 





Table 2. The Parameters of the Cost Function and Generators Limits of the 10 -Unit System 


Unit Prax Pin c b a tup (h) tdown al Shr (h) Scr (h) teold ( init.st 
(MW) (MW) 
1 55 10 1000 40.54 0.12951 8 8 4500 9000 5 8 
2 80 20 950 39.58 0.10908 8 8 5000 10,000 5 8 
3 120 47 900 36.51 0.12511 5 > 550 1100 4 5 
4 130 20 800 39.51 0.12111 5 5 560 1120 4 5 
5 160 50 756 38.53 0.15247 6 6 900 1800 4 -6 
6 240 70 451 46.15 0.10587 3 3 170 340 2 3 
7 300 60 1243 38.30 0.03546 3 3 260 520 2 -3 
8 340 70 1049 40.39 0.02803 1 1 30 60 0 -1 
9 470 135 1658 36.32 0.02111 1 1 30 60 0 -1 
10 470 150 1356 38.27 0.01799 1 1 30 60 0 -1 


Table 3. Load Data for 10 Units, 24h 
Hour Load (MW) Hour Load (MW) Hour Load (MW) Hour Load (MW) 








1 700 ‘it 1150 13 1400 19 1200 
2 750 8 1200 14 1300 20 1400 
3 850 9 1300 15 1200 21 1300 
4 950 10 1400 16 1050 22 1100 
5 1000 11 1450 17 1000 23 900 
6 1000 12 1500 18 1100 24 800 





Table 4. Results of the Economic Dispatching of 10-Unit System 











CM PSO ESCSA 
Pl (MW) 47.8749 47.9492 47.9392 
P2 (MW) 62,1299 62,1460 61,1460 
P3 (MW) 67.8696 67,8343 67,8343 
P4 (MW) 58.8633 58.8639 58.9039 
P5 (MW) 50 50 52 
P6 (MW) 70 70 71 
P7 (MW) 206.8056 206.7936 206.7936 
P8 (MW) 224.5906 224.6132 224.7032 
P9 (MW) 377.5342 377.4816 377.4816 

P10 (MW) 383.4761 383.4313 383.4313 
PL (MW) 49.1165 49.1148 49.1148 
Fuel cost 81115 81115 81113.35 


($/h) 
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Figure 1. Convergence characteristic of fuel cost using Binary ES-CSA for 10-units based UC problem 
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Figure 2. Convergence characteristic of fuel cost using Binary ES-CSA for 10-units based economic dispatch 


6. CONCLUSION 
In this paper, we proposed an eagle strategy based crow search algorithm (ES-CSA) to solve unit 


Economic dispatch and unit commitment problem in smart grid system. The practicality of the proposed 
metaheuristics ES-CSA was tested for 10 power generators test case. The gained results were compared to 
existing results based on PSO and CM methods. It was shown that ES-CSA are superior in obtaining a 
combination of power loads that fulfill the problem constraints and minimize the total fuel cost. ES-CSA 
found to be efficient in finding the optimal power generation loads. ES -CSA was capable of handling the 
non-linearity of ED problem And UC. The evolved power using ES-CSA minimized both the cost of 
generated power, the total power loss in the transmission and maximizes the reliability of the power provided 
to the customers. The programs were developed using MATLAB. The results have shown that our ES-CSA 
to give us a better performance with optimal results in all cases and respecting the constraints imposed. 
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